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Abstract. This paper discusses the properties of the inner puffed-up rim which forms in circumstellar disks when 
dust evaporates. We argue that the rim shape is controlled by a fundamental property of circumstellar disks, 
namely their very large vertical density gradient, through the dependence of grain evaporation temperature on 
gas density. As a result, the bright side of the rim is curved, rather than vertical, as expected when a constant 
evaporation temperature is assumed. We have computed a number of rim models, which take into account this 
effect in a self-consistent way. The results show that the curved rim (as the vertical rim) emits most of its radiation 
in the near and mid-IR, and provides a simple explanation to the observed values of the near-IR excess (the “3 
^m bump” of Herbig Ae stars). Contrary to the vertical rim, for curved rims the near-IR excess does not depend 
much on the inclination, being maximum for face-on objects. We have then computed synthetic images of the 
curved rim seen under different inclinations; face-on rims are seen as bright, centrally symmetric rings on the sky; 
increasing the inclination, the rim takes an elliptical shape, with one side brighter than the other. 


1. Introduction 

The structure of the inner regions of circumstellar disks 
associated with pre-main sequence stars is the subject of 
intense research. Interferometers working in the near in¬ 
frared are providing the hrst direct information on the 
morphology of disks on scales of fractions of AU. They 
show that in the majority of cases the observed visibil¬ 
ity curves are not well reproduced by flared disk models; 
rather, they are consistent with the emission of a ring of 
uniform brightness, of radius similar to the dust evapo¬ 
ration distance from the star (Millan-Gabet et al. I2()()ll 
Tuthill et al. 1200111 . 

The interferometric results provide strong support to 
the idea that the inner disk structure deviates substan¬ 
tially from that of a flared disk because dust evaporation 
introduces a strong discontinuity in the opacity, which 
results in a “puffed-up” rim at the dust destruction ra¬ 
dius, where dust is exposed directly to the heating stel¬ 
lar radiation. The idea of a puffed-up inner rim was pro¬ 
posed by Natta et al. ii^nnTii and developed further by 
Dullemond et al. II2()01I hereafter DDNOl) for Herbig Ae 
stars, to account for the shape of the near-infrared excess 
of these stars (the “3-/im bump”). These authors pointed 
out that the rim had also the right properties to explain 
the early interferometric results of Millan-Gabet et al. 
Il2()()lll . Recent theoretical work by Muzerolle et al. II2()()4II 
has shown that the condition required to produce a puffed- 
up inner rim are indeed likely to exist in most Herbig Ae 
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and T Tauri stars. The concept of such an inner rim has 
been widely used to interpret near-IR interferometric data 
for Herbig and T Tauri stars (Eisner et a, 1. 120041 Muzerolle 
et al. 120031 Colavita et al. 120031 Eisner et al. 120031 Monnier 
and Millan-Gabet |2nn2l Millan-Gabet et al. 12001 1> . Its ef¬ 
fects on the disk structure and emission at larger radii 
have been discussed by Dullemond and Dominik 
who propose that the classihcation of Herbig Ae stars in 
two groups, based on the shape of the far-infrared ex¬ 
cess (Meeus et al. l2001ll . can be interpreted as differences 
between objects where the outer disk emerges from the 
shadow of the inner rim and objects where this does never 
happen. 

In spite of its success in accounting for a variety of 
observations, the actual structure of the rim has not been 
much discussed. DDNOl adopted for their models a very 
crude approximation, namely that the illuminated side of 
the rim is “vertical”, and that its photospheric height is 
controlled by radial heat diffusion behind the rim. Such 
a model, taken at its face value, has the obvious disad¬ 
vantage that the rim emission vanishes for objects seen 
face-on, for which the projection on the line of sight of the 
rim surface is null, and for objects seen edge-on, where 
the rim obscures its own emission. This is clearly incon¬ 
sistent with observations of the SED, which show that all 
the Herbig Ae stars with disks have similar near-IR excess, 
regardless of their inferred inclination (Natta et al. 2001; 
Dominik et al. I2003II . 

The vertical shape of the illuminated face of the rim 
is clearly not physical, as pointed out already by DDNOl. 
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Several effects are likely to “bend” the rim: among them, 
one can expect that radiation pressure on dust grains or 
dynamical instability, due to self-shadowing effects, could 
modify the illuminated face of the rim (Dullemond 2000; 
DDNOl). None of these suggestions, however, has been 
explored further. 

In this paper, we will discuss in detail a different pro¬ 
cess, not mentioned so far, which depends exclusively on 
the basic physics of dust evaporation, i.e., on the de¬ 
pendence of the evaporation temperature on gas density. 
Circumstellar disks are characterized by a very large vari¬ 
ation of the density in the vertical direction, so that the 
dust evaporation temperature varies by several hundred 
degrees in few scale heights; moving vertically away from 
the disk midplane along the rim, dust will evaporate at 
lower and lower temperatures, i.e., further away from the 
central star. This very simple effect curves significantly the 
inner face of the rim, as we will describe in the following. 

The paper is organized as follows. In §2 we describe 
the model we use to compute the rim shape and its obser¬ 
vational properties. The results are presented in §3, where 
we discuss also how the rim depends on dust properties. 
A discussion of the results, in view of the existing obser¬ 
vations, follows in §4, and a summary in §5. 

2. Model description 

Our model of the inner rim of passive-irradiated flaring 
disks joins the two different analytical methods to solve 
the structure of a circumstellar disk proposed respectively 
by Calvet et al. (119911119921 hereafter C92) and Chiang 
and Goldreich CMZl hereafter CG97). The temperature 
in the rim atmosphere is determined using the analytical 
solution of the problem of the radiation transfer as in G92, 
neglecting the heating term due to the mass accretion. 
The vertical structure of the rim is then computed in a 
way derived from GG97 and DDNOl, adding a relation 
between the dust vaporization temperature and the gas 
density as proposed in Pollack et al. (Cnni. As a result 
we obtain a curved model for the inner rim whose features 
are described in the next paragraph. 

Although the expressions for the dust temperature de¬ 
rive from a first order solution of the radiation transfer 
equation, we found (see Appendix) good agreement with 
the correct numerical result in most cases. To zero order, 
we can compute the rim structure avoiding proper radia¬ 
tion transfer calculations. 

In the limit where the incident angle a of the stellar 
radiation onto the disk surface is a <C 1, the equations of 
the temperature for = 0 and 3> 1 (where is the 
optical depth for the emitted radiation) are formally equal 
to the optically thin Tg and midplane Ti temperatures in¬ 
troduced by GG97 (two-layer approximation). Therefore, 
while the two layer approximation is useful to study the 
structure and the emission features (e.g. silicate features 
at lO/rm) of the flaring part of the disk (as in DDNOl), 
it must be abandoned in modeling the inner rim, since 
a ~ 1. Nevertheless we can adapt the GG97 treatment of 


the vertical structure of the disk to the inner rim using 
the appropriate expressions for the dust temperature. 

In order to clarify this concept and to introduce the re¬ 
lation between the vaporization temperature of dust and 
the gas density, the basic equations are briefly summa¬ 
rized. We refer to the cited works for a physical discussion 
of the equations. 

We suppose that the disk is heated only by the stel¬ 
lar radiation and we call a the incident angle beetwen the 
radiation and the disk surface. The incident beam is ab¬ 
sorbed exponentially as it penetrates the dusts and if 
is the optical depth for the emitted radiation, the dust 
temperature T{Td) is given by the relation (G92) 


THrd) = 



fi(2 + + 





( 1 ) 


where /i = sin a and e is an efficiency factor that charac¬ 
terizes the grain opacity, defined as the ratio of the Planck 
mean opacity of the grains at the local temperature T{Td) 
and at the stellar temperature T*, respectively. In writing 
the previous equation we have assumed that the scattering 
of the dust grain is negligible and that the Planck mean 
opacity is defined as: 


KpiT) 


/“ B,{T)Kdv 

S^B,{T)dv • 


( 2 ) 


Following Muzerolle et al. (1200311 . the continuum emission 
of the disk is assumed to originate from the surface char¬ 
acterized by the optical depth = 2/3. 

In the flaring part of the disk, for which /r ^ 1, Eq. o 
can be rewritten in terms of the two layer approxima¬ 
tion, proposed in CG97, in which the interior of the disk 
(with Td^ 1) is heated to the temperature Ti by half of 
the stellar flux, while the other half of the stellar flux is 
reemitted by the superficial layer heated at the optically 
thin temperature Tg'. 


rVd » 1) ^ T/= Tt (3) 

T\Td = 0)^Tt^^-[^^ Tt (4) 


Following DDNOl we assume that the dust component of 
the disk is truncated on the inside by dust evaporation, 
forming an inner hole of radius Rin ■ Inside this hole, only 
gas may exist but as long as the gas is optically thin to 
the stellar radiation (Muzerolle et al. 1200411 we can neglect 
its absorbing effect. Near the inner radius Rin, since the 
rim surface is nearly perpendicularly exposed to the stellar 
radiation (^ ~ 1), Eq. |(2J) and 0 must be replaced by the 
more general relations 

T\Td:>l) = Tt=K‘i + iiie)(^\ Tt, (5) 
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Fig. 1. Width of the region over which 99% of the stellar 
radiation is absorbed as function of e. 


= 0 ) = 



( 6 ) 


derived as limiting cases from Eq. O- 

It is useful to note here that, since the dust in the rim 
is by definition close to the evaporation temperature, the 
value of e is fixed and depends only on the grain absorbing 
cross section. Moreover, to first order, since the difference 
between T(rd = 0) and T{Td = oo) is never very large, 
at the inner radius {fj, = 1) the ratio between these two 
temperature depends only on e and is given by: 


To 


2e + 1 
_ 6(2 + 3e) _ 


1/4 


(7) 


The rim is defined by the condition that the incom¬ 
ing stellar radiation is entirely absorbed by the interven¬ 
ing dust and in the following we will refer to the sur¬ 
face where Ts = 1 to define its shape and location. Once 
the grain evaporation temperature T^yp is fixed, one can 
see from Eq. dU that the rim has a sharp surface when 
6 < Ccr = 1/ VS- In this case the optically thin temper¬ 
ature To of the grains is higher than Too. As soon as a 
grain can survive, the optical depth for the stellar radi¬ 
ation Ts increases rapidly from zero to very high values 
and the rim surfaces characterized by different values of 
the optical depth Tg are all compressed in a vary narrow 
zone, located at Rm, i.e., where Tq = Te„p. For e> e^r, the 
dust temperature increases with the optical depth and the 
transition to optically thick regimes is controlled by the 
geometrical dilution of the stellar radiation; in this case 
the transition from Ts = 0 to very high values occurs in 
a relatively broader region (see Fig. ^|; however, also in 
these cases the rim is located by the Ts = 1 surface with 
an accuracy better than 10-20%. 

With the assumption that the disk is in hydrostatic 
equilibrium in the gravitational field of the central star 


and that is isothermal in the vertical direction z, the gas 
density distribution is expressed by the Gaussian relation 

Pg{r, z) = Pgfi{r) exp(-zV 2 h^), ( 8 ) 


where 



defines the relation between the pressure scale h and the 
interior temperature T^o at a distance r from the central 
star. The temperature Tc is a measure of the gravitational 
field of the central star, expressed by 


Tr = 


GMi,Pg 

kRi, 


( 10 ) 


where pg is the mean molecular weight of the gas. 

Note that the interior temperature of the disk de¬ 
pend on the incident angle a between the stellar radia¬ 
tion and the disk surface, through Eq. ©• The quantity 
p = sin a accounts for the projection of a disk annulus 
on a plane perpendicular to the incident radiation and is 
given by the relation 


f 

r 


dz(r) 

arctan 

A't). 

arctan 

V 

dr 



( 11 ) 


where z{r) is the axisymmetric equation of the surface 
of the disk that we want to determine. For the standard 
flaring model, the incident angle a is given by the relation 


sm a = p~ 


0.4i?* 


d{H/r) 

dr 


( 12 ) 


where H is the photospheric height of the disk, defined 
as the height to which the optical depth of the disk to 
the stellar radiation is Tg = 1, on a radially directed ray. 
For the inner rim we can retain this definition for H and 
place z{r) = H{r) in Ea. itTTll . where H is thus given by 
the relation 


Kp{T.) r ( PI ^ 

- / Pd{z,r)dz = l. 

M JH{r) 


(13) 


As shown in DDNOl, the ratio y = iT/h, between the 
photospheric and pressure height, is a dimensionless num¬ 
ber of the order of 4-6, depending weakly on the angle a, 
the density pd and the Planck mean opacity of dust at the 
stellar temperature T*. Finally, to obtain the dust den¬ 
sity distribution from Eq. © we assume a constant ratio 
Pd/Pg = 0.01. 

We can follow CG97 and DDNOl to obtain a self- 
consistent solution of Eq. (0, 0, ®, ®, 0, CU, C31) to 
determine the structure of the disk, as long as dust evapo¬ 
ration can be neglected. When dust evaporation is impor¬ 
tant, DDNOl have developed an approximate solution of 
the rim/disk structure under the assumption of constant 
Tgyp. The inner rim has a vertical surface toward the star 
located at the dust evaporation distance Rm- The verti¬ 
cal photospherical height H depends on the dust density 
behind the rim; since the rim is higher than the flaring 
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disk, it casts a shadows over the disk. In this shadowed 
region, assuming that there is no external heating except 
the stellar radiation, the pressure height h depends only 
on the radial heating diffusion. The exact determination 
of the structure of the rim in this diffusive region would 
thus require to solve the problem of radiation transport 
in two dimensions (see nnhemond I2()02|> . Since this goes 
well beyond our aims, we adopt the approximated relation 
used in DDNOl 

d{rT^) 

dr h 

Since h oc (from Eq.Q), the temperature Too 

can be eliminated and for h R we obtain the relation 
for h behind the rim 


d{h/r) 1 

dr 8R 


(15) 


The solution of Eq. can therefore be used to determine 
the density Pd(r, z) behind the rim through Eq.®. 

We now introduce the relation between the evapora¬ 
tion temperature of the dust and the gas density. The 
physical reasons for this effect can be easily understood 
thinking of evaporation as the process by which equilib¬ 
rium between the gas pressure and the surface tension of 
the dust grains is reached: the higher the gas density, the 
higher will be the evaporation temperature. We adopt for 
the dust in the disk the model proposed by Pollack et al. 
JTTIMli . In this model, the grains with the higher evapo¬ 
ration temperature are the silicates, which will therefore 
determine the location of the rim. Their evaporation tem¬ 
perature (see Table 3 in Pollack et al.) varies with the gas 
density roughly as a power law, of the kind: 


Tevp=Gp'l{r,Zevp) (16) 

where G = 2000 and 7 = 1.95 • 10“^. Since pg varies expo¬ 
nentially with z (see Eq.®), for Zevp = h the evaporation 
temperature of dust is only 1 % smaller than that on the 
midplane, but for Zgyp = H ^ bh (for a reasonable value 
of x) the difference is about 20 % and dust evaporation 
takes place at a larger distance from the star, curving the 
rim. 

To obtain a self-consistent determination of the struc¬ 
ture of the curved rim, we implemented a numerical 
method able to solve Eq. oni together with the set of equa¬ 
tions @, ®, ®, ®, ®, (HU, (HU. The distance of dust 
evaporation in the midplane is computed for z = 0 in 
the set of equations and is taken as the starting point of 
the radial grid on which the rim structure is computed. As 
a result of the calculations, we obtain the location in the 
(i?, z) plane of the rim surfaces characterized by a con¬ 
stant value of the optical depth. The surface for = 0 is 
thus the evaporation surface of dust grains, for Ts = 1 we 
obtain the surface relative to the photospheric height H, 
defined through the Eq. HU, while for = e • Tg = 2/3 we 
obtain the emitting surface of the rim. 

Moving on the rim surface away from the star, the in¬ 
cident angle a decreases and, when it approaches zero. 


the determination of the rim shape becomes very difficult. 
This is mainly due to the fact that the described solution 
for the radiation transfer neglects the heat diffusion be¬ 
tween contiguous annulus of the rim. Therefore both the 
midplane temperature Too and the pressure height h of 
the rim goes unphysically to zero for a = 0 , according 
to Eq. ® and ® . To avoid this unrealistic behaviour, we 
use the approximated relation discussed previously (see 
Ea. (ll5ll l to determine the pressure height h in the region 
where the diffusion is the dominant heating source. The 
transition distance between the region of the rim heated 
by the star and those heated by the diffusion is determined 
imposing continuity of dh/dr. The photospherical height 
H is then determined as for the vertical rim. 

3. Results 

Using the model described in the previous section, we com¬ 
pute the structure of the inner rim for a disk heated by a 
star with temperature T* = lOOOOK, mass M* = 2.5 Mq 
and luminosity L* = 47 Lq. We take a disk surface density 
E(r) = 2 • 10^{R/AU)~^-^gcm~'^, and a dust-to-gas mass 
ratio dust/gas = 0.01. In our models, this value of E cor¬ 
responds to a midplane gas density of about 10 “®g/cTO^ 
at 0.5AU from the star. Note that the results are not very 
sensitive to the exact value of E, as long as the inner disk 
remains very optically thick. 

The dust properties are those of the astronomical sil¬ 
icates of Weingartner and Draine II2001II . In our models, 
we consider that all the grains have the same size and 
characterize their properties with the quantity e, the ra¬ 
tio of the mean Planck opacity at the evaporation tem¬ 
perature to that at the stellar temperature (see Eq.®). 
The evaporation temperature of silicate grains vary from 
1600K, for gas density of pg = 10~^g/crrA, to lOOOK 
for Pg = (see Ea. dlfill b For the Weingartner 

silicates and a vaporization temperature Te„p ~ 1400 K, 
e ~ 0.08 for grains of radius a = OApm, and grows to 
values of about unity for grain radii > 5pm. 

3.1. The Rim Shape 

The shape of the rim is shown in Fig. lU which plots the 
locus of Ts = 1 (i.e., the photospheric height H) as function 
of R for different values of e. The ending point of the rim 
{RoutI Hout) is when the rim becomes optically thin at the 
stellar radiation. 

For e > Ccr ~ 0.58, corresponding to silicate grains 
bigger than 1.3/rm, the inner radius i?i„, the outer radius 
Rout and the maximum photospheric height of the rim 
Hout, all vary very little, with values ~ 0.50 AU, 
Rout — 0.65 AU and Hout — 0.09 AU. For smaller values of 
e, the rim gets steeper and the inner radius increases. For 
e = 0.08, corresponding to grains with radius 0.1pm, the 
rim has Rin = 1-16 AU, Rout = 1-34 AU and Hout = 0.14 
AU. 

The shape of the rim can be roughly characterized by 
the ratio of its maximum height Hout over width Ai? = 
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Fig. 2. Photospherical height H (ts = 1) of the inner rim for different values of e, as labelled. The curves have been 
computed for a star with T* = lOOOOK, = 2.5Mq, L* = 47Lq. Each curve ends where the rim becomes optically 
thin to the stellar radiation. 



Fig. 3. Effective temperature (i.e., T{Td = 2/3)) along the 
rim surface for different values of e, as labelled. Each curve 
is plotted as function of R/Rm, where Rm is the distance 
of the rim from the star on the disk midplane. Note that, 
in fact, not only R but also z changes along each curve, 
as shown in Fig. |21 For e > €„ — 0.58 T^ff is almost 
independent of the grain opacity. 

{Rout — Rin)- This quantity, which is nominally infinity in 
a vertical rim, becomes in the curved rim models ~ 0.6 
for e > Ecr and is ^ 0.8 for e ^ 0.14 — 0.08. In other words, 
as grains grow, the inner rim approaches the star but the 
bending of the surface varies very little. 

For any given value of e, the dust temperature along 
the rim is not constant (as in the vertical rim) but de¬ 
creases from values of about 1400 K, typical of silicate 
evaporation temperatures at density of 10~^g/cm^, to 
1200 K for density of 10~^^5/cm^. Fig. |21 plots the ef¬ 
fective temperature of the rim (i.e., T{Td = 2/3)) along 
its surface for different values of e. For e > Ccr the effective 
temperature is equal to the vaporization temperature; for 
e < Ccr the effective temperature is generally lower since 
the dust temperature decreases with the optical depth, ac¬ 
cording to Eq. O- This variation, however, is sufficiently 
small that most of the rim emission occurs in the near-IR, 
as for the vertical rim. 


3.2. The Rim SED 

The fraction of stellar luminosity intercepted by the rim, 
given by the ratio Hout/Rout, varies from 10%, for small 
values of e, to 14%, for e > €cr- Assuming that the rim is 
in thermal equilibrium, the intercepted radiation is equal 
to the total emitted flux For each e, we have com¬ 
puted the spectral energy distribution of the rim emission 
for different inclination angles. As discussed in the rim 
emission is computed as that of a blackbody at the local 
temperature along the rim = 2/3 surface. We will come 
back in the Appendix to this assumption, which, in any 
case, gives a very good approximation to the global prop¬ 
erties of the rim emission (see also Muzerolle et al. 120031 
C92). Most of the emission, as expected from the range of 
temperatures, occurs in the near-IR. We have computed 
the fraction of the stellar luminosity re-emitted by the rim 
in the wavelength range 1.25-7.0 gm for different values 
of e as function of the inclination (Fig. 0)1. This near-IR 
excess peaks at zero inclination, where has values between 
~ 10% (smaller grains) and ~ 20% (larger grains). As the 
inclination increases, the near-IR excess decreases slowly, 
reaching values between 5% and 8%, depending on e. For 
inclination higher then ~ 80° the rim emission is self- 
absorbed. Note that for the large grains, with e > ecr, the 
near-IR excess becomes almost independent of e. 

The behaviour described in Fig. Elis very different from 
what one obtains in the case of the vertical rim. As de¬ 
scribed in the previous section, neglecting the dependence 
of the evaporation temperature of grains from gas density 
results in the inner face of the rim to be vertical (as in 
DDNOl). With e = ter and T^vp — IdOOAT (equal to the 
evaporation temperature at the inner radius of the curved 
rim), the inner radius Rin is the same as for the curved 
rim; the photospheric height, evaluated using the approxi¬ 
mation described by Eq. (Cnj, is also the same {Hout =0.09 
AU); the fraction of stellar luminosity intercepted by the 
vertical rim is 17%. However, the value of observed 

for different inclination angles is very different, with a very 
strong (and opposite) dependence on i. (Fig. EJ. In par¬ 
ticular, Finin vanishes for face-on disks, and is maximum 
(of the order of 20%) for very inclined systems, just before 
self-absorption sets in. 
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There are also differences between the curved and ver¬ 
tical rim models concerning the shape of the predicted 
3/xm bump. While the vertical rim has a constant tem¬ 
perature of 1400K over all its surface, in the curved rim 
the temperature varies from 1400K on the midplane to 
about 1200K at the outer edge (see Fig. 121 and the SED 
is broader than a single-temperature black body. In prac¬ 
tice, however, this is only a minor effect. 

3.3. Rim Images 

Fig- El shows synthetic images of the curved rim at 2.2/rm 
for grains with e = 0.08 (radius of about 0.1 /im) and e = 
Ccr = 0.58 (radius of about 1.3 ^m). As described in §3.1 
(see Fig. Inland Fig. 121), the inner radius of the rim is larger 
for smaller grains and the surface brightness is lower, due 
to the lower effective temperature of the emitting surface. 

For the comparison with the vertical model of the rim, 
the left panels of Fig. |21 show the images of the vertical 
rim, calculated for e = 0.58. The largest difference is at 
low inclinations: for i = 0° (face-on disks) the vertical rim 
vanishes, as the projected emitting surface along the line 
of sight is zero, while the curved rim has a centrally sym¬ 
metric ring shape. For a distance of 144pc, the inner radius 
of the ring, for e = ecr, is about 4 milliarcsecond (mas); 
its brightness peaks practically at Rin (Rpeak/Rin ~ 1.05) 
and decreases slowly (by about a factor of two) outward, 
until at Rout it drops to very low values. At this distance, 
the width of the rim (roughly its FWHM) is about 0.8 
mas. 

For higher inclination the central symmetry is lost and 
the projected image is an ellipse with one edge brighter 
than the other. In general, however, the brightness distri¬ 
bution of the rounded rim is much more symmetric than 
that of the vertical rim. 

3.4. Grain size distribution 

The results shown so far have been computed assuming 
that all the grains have the same composition and size. 
This is in practice unrealistic, and one wonders how the 
results will change if grains with different properties are 
present. Different grains absorb differently the stellar ra¬ 
diation and reach very different equilibrium temperatures. 
The degree of complexity of the radiation transfer prob¬ 
lem increases remarkably (see Wolfire and Cassinelli ll 9861 
nm and there is no approximate solution such as Eq. 0 ). 
In fact, to the best of our knowledge even numerical solu¬ 
tions are not available to describe accurately the transition 
region where some grains are cooler than evaporation and 
survive while others do not. 

We can however estimate the effect of a grain mixture 
in the following way. Let us consider for simplicity a classi¬ 
cal MRN grain size distribution, characterized by a power 
law a“'^ with a the radius of the dust grains, <7=3.5, and 
a varying between a minimum value Umm = 0.01 gm and 
a maximum value Umax ^ Umin- From what we have dis- 



Fig. 4. Near infrared emission of the curved rim (inte¬ 
grated between 1.25 — 7fim and normalized to L*) for dif¬ 
ferent inclination angles and different values of e. Starting 
from the bottom, the curves refer to e = 0.08 (long- 
dashed), 0.14 (short-dashed), 0.25 (dotted), 0.58 = Ccr 
(short-dash-dotted), 0.92 (long-dash-dotted). For the last 
two values of e the emission is almost the same and the 
lines overlap. The continuum line plots the emission of the 
vertical rim, computed for e = Ccr- 

cussed in (121 if Umax is smaller than 1.3gm (corresponding 
to e = Ccr for the adopted silicate grains), all the grains 
have a maximum temperature at = 0 and the higher 
is the value of e, the smaller is the evaporation distance 
from the star. As soon as the largest grains can survive, 
the stellar radiation is rapidly absorbed and all the other 
grains will also survive. In this case, the shape of the inner 
rim is controlled by the largest grains in the distribution. 
The shapes in Fig. Hand the emitted fluxes shown in Fig. 
E] should not change significantly as long as one interprets 
e as relative to the largest grains. 

The situation is more complex if Umax > t.Sfim, since 
the grains with e > Ccr can survive near the star only in 
an optically thin regime. However, our single-grain models 
show that there is very little difference in the location and 
shape of the rim as soon as e > Ccr- This suggests that 
varying Umax above the “transition” value (i.e., the value 
of e at which the dependence of T on r changes from de¬ 
creasing to increasing) will not change the rim properties 
any further, at least to zero order. 

It is also likely that grains have not just different size 
but also different chemical composition. In this case, the 
rim location and properties are determined by the dust 
species with the highest evaporation temperature, as long 
as its contribution to the disk opacity is sufficient to make 
the disk optically thick. As discussed in §2, in the Pollack 
et al. (1994) dust model silicates have the highest evapo¬ 
ration temperature. This is why we have performed all our 
calculations for silicate grains. However, one should keep 
in mind that, in some dust models, graphite contributes 
most of the opacity at short wavelengths; since graphite 
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Curved Rim - E = 0.08 


Curved Rim - E = 0.58 - Vertical Rim 



Distance from the central star (milli arcsec) 


Fig. 5. Synthetic images of the curved rim for different values of e and inclination angle (z = 0° is face-on); comparison 
with the vertical rim. The left and central panels show the images of the curved rim calculated for e = 0.08 (silicate 
grains with radius a = 0.1/rm) and e = Ccr = 0.58 (a = I.3/rm). The right panels show the images for the verical rim 
calculated for e = 0.58. The surface brightness of the rim, plotted in colors, is computed for a wavelength of 2.2/iTO. 
The stellar parameters are as in §3, the distance is d = 144pc. The vertical rim for z = 0° has zero brightness. Note 


that as discussed in §3 for e > ecr the images of the curved 

has a much higher evaporation temperature {Tf.yp ~ 2000 
K) than silicates, its presence would move the rim much 
closer to the central star. Although the details of the rim 
shape may change, its curving, which is caused only by 
the dependence of Tevp on the density, will not disappear. 

4. Discussion 

4.1. The 3fj,m bump 

The curvature of the inner side of the rim has important 
consequences on the observable near-IR excess, as shown 
in Fig. ^ which shows that the predicted near-IR excess 
ranges between about 10 and 20%, and that there is no 


rim remain the same. 

significant dependence on the inclination of the disk with 
respect to the observer. These results compare very well 
with the existing observations. Fig. M shows the observed 
values for a well-studied sample of Herbig Ae stars from 
Natta et al. (2001) and Dominik et al. (2003); of a total 
of 16 objects, only one (HD 142527) has LMin/Li, > 0.25, 
and one (HD 169142) < 0.09. 

The observations of the near-IR SEDs of Herbig Ae 
stars have also shown that there is no systematic varia¬ 
tion of the near-IR excess with the inclination of the disk. 
The sample of Fig. El includes some rather face-on ob¬ 
jects such as AB Aur (Ljv/fl/A* ~ 0.20, i ~ 20° — 40°; 
Fukagawa et al. IWill and HD 163296 {Lnir/L^ ~ 0.21, 
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Fig. 6. Histogram of the observed values of LMm/Li, for 
a sample of 16 Herbig Ae stars from Natta et al. (2001) 
and Dominik et al. (2003). Note that the Dominik et al. 
values are computed between 2 and 7/rm, while those from 
Natta et al. between 1.25 and 7 /rm. The uncertainties are 
in most cases rather large, due to the subtraction of the 
photospheric flux, which contributes significantly to the 
observed ones at short wavelengths, and to variability (see 
Natta et al. for some examples). 


i ~ 30°; Grady et al. r2000|l and some UXORS variable, 
such as UX Ori, WW Vul and CQ Tau, which are gen¬ 
erally considered to be close to edge-on, and have values 
of LMiR/Li, ~ 0.12 — 0.25. This, which has been a puzzle 
for the vertical rim, finds a natural physical explanation 
in our models, which predict that the near-IR excess de¬ 
pends little on the inclination and that, in particular, does 
not vanish for face-on objects. Actually, our models pre¬ 
dict that the largest excesses should be seen in face-on 
objects, and it would be interesting to explore in detail if 
this is indeed the case. However, for such a study to be 
significant, one would need a large sample of objects with 
known inclinations, which is at present not available. 

Another interesting aspect of our results is the depen¬ 
dence of on grain properties. If taken at face- 

value, one should expect that only objects with relatively 
large grains (greater than about I /im) can have large val¬ 
ues of Lnjr /and that low values of the near-IR excess 
can only occur in face-on objects with small grains. These 
properties of the rim emission are potentially important 
for a better understanding of grain properties in the inner 
disk and should be pursued further in the future, com¬ 
bining observations of disks with well known inclinations 


with models that explore a broader range of grain prop¬ 
erties than we considered in this paper. 

4.2. The rim radius 

One side product of our models is i?i„, the distance from 
the star to the rim on the disk midplane. For a face-on 
disk, Rin practically coincide with the position of the peak 
of the near-IR brightness. For fixed values of the stellar 
parameters, Rm depends on the dust properties; its value 
is shown in Fig. 0 for different values of e. Larger grains 
can live closer to the star than smaller grains and, be¬ 
cause of the inversion in the temperature gradient, Rin 
does not change much with the grain size once this ex¬ 
ceeds the value for which e = Ccr, so that for L* = HLq, 
Rin ^ 0.5AU. Rin scales roughly as L* so that one can 
expect Rin ^ O.I AU for L* ~ ‘2Lq. These values, espe¬ 
cially for small grain sizes, are larger than the predictions 
of simple, optically thin calculations, because Eq.(I) cor¬ 
rectly includes the effect of the diffuse radiation field. Note 
that if backward scattering cannot be neglected, the tem¬ 
perature at T = 0 will be even higher, increasing further 
the size of the inner hole (see Appendix). 

Values of Rin have been derived for a number of objects 
from near-IR interferometric observations ( see, for exam¬ 
ple, Akeson et al. I2()()()l I2002L Millan-Gabet et al. I200IL 
Monnier and Millan-Gabet Em Colavita et al. Em 
Eisner et al. Em Em Em. The results depend not 
only on the assumed disk model, but also on other quan¬ 
tities such as the stellar SED, the disk inclination etc., 
so that a direct comparison with our values of Rin can 
be misleading. We note, however, that Rin ~ 0.5 AU for 
intermediate-mass objects is roughly consistent with the 
observations, suggesting that in many objects grains have 
grown to sizes a > 1.3/im (see also Monnier and Millan- 

Gabet Em. 

This conclusion, of course, needs to be taken with great 
caution. There are a number of effects that can change the 
model-predicted Rim in addition to different grain sizes. 
For example, on one side smaller values of Rin can be 
due to some low-density gas in the inner disk hole, able 
to absorb the UV continuum from the star (Monnier and 
Millan-GabetEnnH Akeson at al. 1200511 . Also, the presence 
of same graphite, which has an evaporation temperature 
of about 2000 K, can decrease Rin- On the other side, the 
presence of accretion at a significant rate can increase Rin 
by heating grains in the inner disk to temperatures higher 
than those produced by the photospheric radiation alone 
(Muzerolle et al. 120041) . 

In practice, one needs to compare in detail model pre¬ 
dictions of the quantities measured with the interferom¬ 
eters (visibility curves, their dependence on baseline and 
hour angle, phase measurements) for specific, well known 
objects. We are currently performing such a study with 
the aim of assessing the impact of curved rim models on 
the interpretation of current (and future) interferometric 
observations. In particular, we will include a discussion 
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of the effect of the asymmetries of the rim projected im¬ 
ages when not face-on. The results will be presented in a 
forthcoming paper (Isella et al. 2005). 

4.3. LkHalOl 

To the best of our knowledge, the only image of the inner 
region of a circumstellar disk is that obtained by Tuthill et 
al. (EDniJ, for LkHalOl using Keck in the H and K bands. 
The results are reminescent of our images, showing an 
elliptical ring with a side much brighter than the other. 
The similarity is very interesting but one should keep in 
mind that LkHalOl is an early B star with a luminosity 
between 500 and 50000 Lq, depending on its distance. Our 
model may not apply to such an object (e.g., Monnier et 
al. 1200511 . 

5. Summary and Conclusions 

In this paper, we discuss the properties of the inner puffed- 
up rim which forms in circumstellar disks when dust evap¬ 
orates. The rim existence has been argued for, starting 
from the work of Natta et al. (2001), both on theoretical 
and observational grounds. Here, we investigate the shape 
of the illuminated face of the rim. We argue that this shape 
is controlled by a fundamental property of circumstellar 
disks, namely their very large vertical density gradient, 
through the dependence of grain evaporation temperature 
on density. As a result, the bright side of the rim is nat¬ 
urally curved, rather than vertical, as expected when a 
constant evaporation temperature is assumed. 

We have computed a number of rim models, which 
take into account this effect in a self-consistent way. A 
number of approximations have been necessary to perform 
the calculations, and we discuss their validity. We think 
that the basic result (i.e., the curved shape of the rim 
illuminated face) is in fact quite robust. 

For a given star, the rim properties depend mostly on 
the properties of the grains, and very little on those of the 
disk itself, for example the exact value of the surface den¬ 
sity. The distance of the rim from the star is determined 
by the evaporation temperature (at the density of the disk 
midplane) of the dust species that has the highest evap¬ 
oration temperature, as long as its opacity is sufficient to 
make the disk very optically thick; in the model of Pollack 
et al. of the dust in accretion disks, silicates have 

the highest evaporation temperature. Therefore, we have 
assumed in our models dust made of astronomical sili¬ 
cates, and varied their size over a large range of values. 
We find that the rim properties do not depend on size as 
soon as a>1.3/rm; the values of the rim radii observed with 
interferometers suggest that in many pre-main-sequence 
disks grains have grown to sizes of 1-few /rm at least. 

The curved rim (as the vertical rim) emits most of its 
radiation in the near and mid-IR, and provides a simple 
explanation to the observed values of the near-IR excess 
(the “3 fj,m bump” of Herbig Ae stars). Contrary to the 
vertical rim, for curved rims the near-IR excess does not 



e 


Fig. 7. Behaviour of the inner radius of the rim for differ¬ 
ent values of e, evaluated for the star and disk parameters 
described in m Note that for e > ecr ~ 0.58 (correspond¬ 
ing to silicate grains of 1.3/rm) the radius of the inner rim 
is almost constant around 0.5AU. 


depend much on the inclination, being maximum for face- 
on objects and only somewhat smaller for highly inclined 
ones. This is in agreement with the apparent similarity of 
the observed near-IR SED between objects seen face-on 
and close to edge-on. 

Finally, we have computed synthetic images of the 
curved rim seen under different inclinations. Face-on rims 
are seen as bright, centrally symmetric rings on the sky; in¬ 
creasing the inclination, the rim takes an elliptical shape, 
with one side brighter than the other. However, the bright¬ 
ness distribution of curved rims remains at any inclination 
much more centrally symmetric than that of vertical ones. 
In a forthcoming paper (Isella et al. 2005) we will discuss 
the application of the curved rim models to the interpre¬ 
tation of near-IR interferometric observations of disks. 
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Appendix A: Appendix 

In this appendix we discuss in detail some of the assump¬ 
tions on radiation transfer made in building our rim mod¬ 
els. We use as templates the results of a radiation transfer 
code developed by E. Kriigel, which is described in Habart 
et al. (2004). The code considers a plane-parallel slab of 
dust, illuminated on one side by a star with properties as 
in §3. 
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Fig. A.l. Temperature as function of the optical depth to 
the stellar radiation Tg for small grains (e = 0.08) and very 
large grains (e = 1). The dashed curves are the results of 
the radiation transfer code, the solid lines the temperature 
obtained from Eq. O- 

A.l. Evaluation of Eq.(l) for single grains 

We show in Fig. ro the comparison of the temperature 
derived from Eq.(l) with the results of the radiation trans¬ 
fer code for two cases, one with e = 0.08 (small grains) and 
one with e = 1 (very large grains). In both cases scatter¬ 
ing is neglected. The temperature is plotted as function 
of the optical depth to the stellar radiation y,. One can 
see that Eq.(l) gives the correct values of T for Ts = 0 
and for Tg = oo. For intermediate values of Tg, the results 
are in both cases accurate within 10%. This is more than 
adequate for the purposes of this paper. 

A.2. SED 

We compare now the SED of a black-body at temperature 
Tigd = 2/3) from Eq.Q, adopted in this paper to describe 
the rim emission (see §3.2) to the results of the radiation 
transfer code (dashed curves) for the same grains of Fig. 
mu The approximation is very good for large grains. For 
small values of e, the discrepancy is larger, and, in par¬ 
ticular, one cannot reproduce any dust feature. However, 
also in the case e = 0.08, the difference in Ljv/_r/L* is only 
of 4%. 

A.3. Scattering 

Scattering has the obvious effect of increasing the grain 
temperature at the slab surface and decreasing it at large 



Fig. A.2. Comparison of the SED computed with the ra¬ 
diation transfer code (dotted lines) and as a black-body 
at T{Td = 2/3) from Eq.(l) (solid line) for e = 0.08 and 
e = 1.0, as labelled. 

optical depth. We compare in Fig. roi the temperature 
profile of grains with e = 0.08 when scattering is prop¬ 
erly included in the radiation transfer (dashed line) rather 
than suppressed (Qsc = 0.). The curved are labelled with 
the values of uj, defined as Qgc{l — g)/Qsc{l — g) + Qabs, 
where g measures the asymmetry of the scattering phase 
function; J; = 0 if Qabs = 0 or if the scattering is forward 
peaked {g = 1). All quantities are averaged over the stellar 
radiation field and the corresponding SEDs show that the 
value of Ljv//j/L* decreases when scattering is included 
by about 25%. 
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